EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



1 July, 2002 



A Study of Inclusive 
Double— Pomeron- Exchange 
in pp — > pXp at y/s = 630 GeV 

A. Brandt 1 , S. Erhan a , A. Kuzucu 2 , M. Medinnis 3 , 
N. Ozdes 2 ' 4 , RE. Schlein b , M.T. Zeyrek 5 , J.G. Zweizig 6 
University of California*, Los Angeles, California 90024, U.S.A. 

J.B. Cheze, J. Zsembery 
Centre d'Etudes Nucleaires-Saclay, 91191 Gif-sur-Yvette, France. 

(UA8 Collaboration) 

Abstract 

We report measurements of the inclusive reaction, pp — » pXp, in events where 
either or both the beam-like final-state baryons were detected in Roman-pot spec- 
trometers and the central system was detected in the UA2 calorimeter. A Double- 
Pomeron- Exchange (DPE) analysis of these data and single diffractive data from the 
same experiment demonstrates that, for central masses of a few GeV, the extracted 
Pomeron-Pomeron total cross section, cr^-p, exhibits an enhancement which ex- 
ceeds factorization expectations by an order-of-magnitude. This may be a signature 
for glueball production. The enhancement is shown to be independent of uncertain- 
ties connected with possible non-universality of the Pomeron flux factor. Based on 
our analysis, we present DPE cross section predictions, for unit (1 mb) Pomeron- 
Pomeron total cross section, at the Tevatron, LHC and the 920 GeV fixed-target 
experiment, HERA-B. 
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1 Introduction 



We study the Double-Pomeron-Exchange (DPE) "diffract ive" process|IJ, depicted in 
Fig. |l](a), using the reaction: 

pp — > p X p. (f ) 

When the final state p and p both have large Feynman-Xp, the process proceeds with the 
exchange of two (virtual) gluon-rich colorless systems called Pomerons. These systems, 
which carry a small fraction of the beam momentum of the two approaching hadrons, 
Ap/p = £ = 1 — x p , collide and constitute the entire effective interaction between the two 
beam particles. This leads to the presence of "rapidity-gaps" , or regions of pseudorapidity 
with no particles between the outgoing p and p and the central system, X. The system 
X with invariant mass, Mx, is the result of the Pomeron-Pomeron interaction. To good 
approximation, Mx 2 = s' = ^i^s (we use the symbols, Mx 2 and s', interchangably) ; 
thus, a given Mx is produced at smaller £ values when the cm. energy is larger. Single 
diffractive processes (see Fig. |l|(c)) appear to be essentially pure Pomeron-exchange 
when £ < 0.03. We expect this to be also the case in DPE. 

The DPE process is the closest we can come to pure gluon interactions. As such, 
it may be a splendid glueball production process ||. At the very high energies of the 
LHC, "diffractive hard scattering" in React. [1] may have advantages as a relatively clean 
production mechanism of rare states. Diffractive hard scattering was proposed in Ref. j|] 
and discovered in pp interactions by the UA8 experiment j5| at the SppS-Collider (y/s = 
630 GeV) and in ep interactions by the HI and ZEUS experiments at HERA|| (see also 
Refs. ^ for studies of these effects at the Tevatron). First studies of hard scattering 
in React. p] were made at the SppS-Collider and at the Tevatron [pi], |TT| . Based on 
UA8's observation of the "super-hard" Pomeron in single diffractive dijet production ||, 
a small fraction (^ 10%) of all hard-scattering DPE events at the LHC may be high-mass 
gluon-gluon collisions. 

In the present paper, we present final results on React. [I] from the UA8 experiment at 
the CERN SppS-Collider. UA8 was the first experiment in which data acquisition from 
a large central detector was "triggered" by the presence of outgoing beam-like protons or 
antiprotons. The final-state baryons were measured in UA8 Roman-pot spectrometers [|12 



which were installed in the outgoing arms of the same interaction region as the UA2 



experiment ||13|| : the central system, X, was measured in the UA2 calorimeter. At 630 
GeV center-of-mass energy, Mx = 6.3 (18.9) GeV when £i and £2 are both equal to 0.01 
(0.03). 

Previous measurements of React. [I] with exclusive final states were made in pp 



interactions at the CERN Intersecting Storage Rings with cm. energy, i/i = 63 GeV, 



and with a beams [15] at \/s = 126 GeV. The advantage of the present ten-times higher 
cm. energy is that much smaller values of £ are accessible for a given produced Mx, 
thereby enhancing the purity of the Pomeron-exchange component. 

As seen in Fig. [I], Reaction [l] is intimately related to the inclusive single-diffractive 
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reactions: 

pp — > pX or pp — > Xp , pp — > X p. (2) 

In Reactions |2j, the Pomeron from one beam particle interacts with the second beam 
particle. Fig. |T|(c) depicts the single diffractive reaction, while Fig. |I](b,d) show the cor- 
responding (dominant) Triple- Regge diagrams of DPE and inclusive single diffraction. 

At small four-momentum transfer^], \t\, triple-Regge predictions of inclusive single 
diffractive cross sections are found to increase more rapidly than do the observed cross 
sections JTB], [17| and violate the unitarity bound above s/s = 2 TeV. The observed large 



damping effects in the data are believed to be due to multiple-Pomeron-exchange effects, 
which phenomenologically are equivalent to a smaller effective Pomeron trajectory inter- 
cept with increasing energy (ff8l |T9. However, despite these unitarizing effects, effective 



vertex factorization appears to remain valid to an astonishing degree fll9[ . In the present 
analysis, we assume its validity. 

In terms of the Triple-Regge model, the cross section for React. |2] may be written 
as the product of the Pomeron-proton total cross section, with the flux factor for 
a Pomeron in the proton , F-p/ p (t,^). Since it is our working assumption that the same 
FT/p(ty describes the Pomeron-proton vertices in both Reacts. || and [l], the cross section 
for React. [I] is given by the product of the Pomeron-Pomeron total cross section, <T^p(s'), 
with two flux factors. See, however, the discussion in Sect. [7] on systematic uncertainties 
due to a possible non-universality of F-p/ p (t,£). The essential result of this paper will be 
shown to be insensitive to such effects. 

The empirical F-p/ p (t,£,) has been "fine-tuned" in fits of the following equation to all 
available data on React. |] at the SppS p| and ISR fl20| : 



= F v/p (t,Q-o%W = [X .|F 1 (t)| a . e ".^*)].[a .((^ M + i2(^-^)]. (3) 



|Fi(t)| 2 is the Donnachie-Landshoff [ 2TJ form factorf]. The right-hand bracket in Eq. |3], the 
Pomeron-proton total cross section, is assumed to have the same form that describes the s- 
dependence of real particle total cross sections. The best values of the fitted parameters [Q] 
in Eq. | aref]: 

Ka = 0.72 ±0.10 mb GeV" 2 

R = 4.0 ±0.6 

b = 1.08 ±0.20 GeV- 2 



The effective Pomeron trajectory is found [|19[] to be s-dependent and, at the energy of 
the UA8 experiment (a/s = 630 GeV), is: 

a(t) = 1 ± e + a't + a"t 2 = 1.035 + 0.165* + 0.059t 2 (4) 



throughout this paper, we refer to the squared four-momentum transferred to the final-state p or p 
as the "momentum transfer" 
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The fits are also consistent with the existing CDF results at the Tevatron |T^, |l^] 



2 



while, over the ISR energy range (s = 549 to 3840 GeV 2 ): 



e(s) = (0.096 ± 0.004) - (0.019 ±0.005) ■ log(s/549). 
a'(s) = (0.215 ±0.011) — (0.031 ± 0.012) • %(s/549). 
a"(s) = (0.064 ± 0.006) - (0.010 ±0.006) • %(s/549). 



The quadratic term|| in a(t) corresponds to a "flattening''^, or departure from linear 
behavior, of the effective Pomeron trajectory at high-|t|. Direct evidence for this flattening 
of the trajectory can be obtained by looking at the behavior of the UA8 single diffractive 
data at large- \t\. Figure || shows the observed Feynman-x p distributions for different 
bands of \t\ between 1 and 2 GeV 2 . Since the geometrical acceptance |2J depends linearly 
and weakly on x p in this figure, the pronounced peaks near x p = 1 reflect the physics of 
diffraction and are seen to persist up to \t\ of 2 GeV 2 . They are due to the (approximate) 
I/My behavior of Triple- Regge phenomenology. If the trajectory did not flatten, but 
continued to drop linearly, the diffractive peak would tend to disappear. For example, 
with a trajectory, a{t) = 1.08 ± 0.25t, the peak would disappear at —t = 2.3 GeV 2 
(corresponding to: 2a (t) — 1 = 0). Thus, the persistence of the diffractive peak in Fig. || 
is the most direct evidence that the effective Pomeron trajectory flattens at large- 1 1 1 . 

The question arises as to whether Pomeron-exchange is still dominant for \t\ > 1 GeV 2 , 
where most of the data in the present experiment are. The self-consistency of our Triple- 
Regge analysis in describing both single-diffraction and double-Pomeron-exchange data 
is one supporting argument. Another important point is that the set of all SppS and 



ISR high-|t| data agree [17] with a "fixed-pole" description without damping. Another 
argument is that the hard Pomeron structure found in the UA8 jet event analysis M 
is consistent with that found in the analysis of low-|t| data at HERA p3| . Thus, our 



working assumption is that Pomeron-exchange dominates React. [I] in the momentum- 
transfer range, 1 < —t < 2 GeV 2 . Based on the results of earlier studies0 of diffraction, 
we can ignore Reggeon exchange when £ < 0.03. 

The differential cross section for the DPE process, React. [I], is: 

d 6 <JDPE 



dC, 1 dC,2dtidt 2 d(f)id(f)2 



F P/p (t 1 ,Z 1 )-F v/p (t 2 ,S 2 )-<T%,(s'). (5) 



The variables, (£i,ti,fa), describe each of the emitted Pomerons at the outer vertices in 
Fig. |l|(a), which are uniquely given by the measurement of the associated outgoing p (or 
p) in the final state. Although there is no explicit 0-dependence on the right- hand- side 
of Eq. |5] and the Pomerons are emitted independently and isotropically, <fi correlations do 
result, because significant regions in the 6-dimensional space, t\, (pi, t 2 , fa), are 
unphysical and give s' < 0. This point is discussed further in Sect. |] in connection with 
Monte-Carlo generation of events according to Eq. [| 

Using Eq. [|, our goal is to extract cr^-p from our data on React. [I] and to determine its 
energy (s') dependence. In particular, we wish to know whether there are enhancements 

4 This flattening is also claimed to be seen by the ZEUS experiment [^2) at DESY in photoproduction 
of low- mass vector mesons (p° and 0°) 
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at small s' which could be due to a strong Pomeron-Pomeron interaction and possible 
glueball production. In the large-s' region where Pomeron-exchange dominates, a™, is 
related by factorization to the Pomeron-proton and proton-proton total cross sections: 

tot I i\ l a Vp( S ')} 2 / n \ 

-vvC) = (6) 

This is seen with reference to the ratios of forward elastic amplitudes for the three pro- 
cesses shown in Fig. ||. A generalized optical theorem |^4|. |25|] then leads to Eq. |6| between 



total cross sections. Ryskin p6| has pointed out that the three cross sections must be 
evaluated at the same value of s'. 

Despite the fact that the cross sections, <Xp* and cr^p, can only be extracted from 
data in product with the constant K in -Fp/ P (t,£) (K 2 in the case of Op V and K in the 
case of dp*; see Eqs. |3] and §), we see m Eq. |6] that such factors of K cancel. Thus, the 
factorization test does not require knowledge of K. However, absolute values of either cr^* 
or Opp can only be given for an assumed value of K, for example by using the Donnachie- 
Landshoff model |?T| with K = 9f3 2 /(4ir 2 ) = 0.74 GeV -2 , which arises from an analysis of 
elastic scattering data. Although different multi-Pomeron-exchange effects in diffraction 
and elastic scattering mean that this value of K is only approximate, we nonetheless do 
quote values for a^j, in the closing sections of this paper, assuming K = 0.74 GeV -2 . 

After describing the experiment and the event selection, we discuss the Monte-Carlo 
event generation of React. [I], and the determination of a^ v . In the process, we com- 
pare the results from two different data samples, one in which both p and p are detected 
and the other in which only the p or the p is detected. We then test the factorization 
relation and come to our conclusions which demonstrate, among other things, an over- 
all self-consistency of our phenomenological description of single diffraction and double- 
Pomeron-exchange. Based on these conclusions, we calculate predictions for double- 
pomeron-exchange yields at the Tevatron and LHC and also at the HERA-B fixed-target 
experiment. 



2 Apparatus & trigger 

Detailed descriptions of the UA8 apparatus, its properties, triggering capabilities and 
interface to the UA2 experiment (TB| are given elsewhere [|T^j. Thus we only provide a 
brief summary of the spectrometer here. 

UA8 constructed Roman-pot spectrometers [|T^] in the same interaction region as UA2, 
in order to measure the outgoing "beam-like" p and/or p in React. |I]or React. 0, together 
with the central system using the UA2 calorimeter system fpT3|| - There were four Roman- 
pot spectrometers (above and below the beam pipe in each arm) which measured p and/or 
p with x p > 0.9 and 0.8 < \t\ < 2.5 GeV 2 . Fig. (| shows one spectrometer with the 
trajectories of 300 GeV particles (x p ~ 0.95) emerging from the center of the intersection 
region with minimum and maximum acceptable angles (solid curves). The lower (upper) 
edge of the shaded area corresponds to the minimum (maximum) accepted angles of x p = 1 
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tracks. The trajectory corresponding to the lower edge of the shaded region is 12 beam 
widths (a) from the center of the circulating beam orbit. 

Particle momenta in the Roman pot spectrometers were calculated in real-time by a 
dedicated special-purpose processor system [[12], ^7|, thereby providing efficient low-rate 
p and p triggers. Improved final-state proton and/or antiproton momenta are calculated 
offline using the reconstructed vertex position (if it exists), given by the UA2 central 
chamber system fI5| , and points reconstructed from hits in Roman pot chambers 1, 2 and 
3. Chamber 4 was also used in the fit, if a track traversed it. 

Figure [5] shows a "beams-eye" view of the UA8 chamber aperture which is closest to 
the center of the interaction region. The four-lobed curve in the figure illustrates the 
contour of the beam pipe which matches that of the quadrupole-magnet pole pieces. The 
overlap between the beam pipe and rectangular chambers above and below the beam 
illustrates the limited azimuthal ranges in which a particle may be detected. These are 
centered at ~ 90° and ~ 270°. Data were recorded with the bottom edge of each pot 
set, in different runs, at either 12 beam widths (12a) or 14cr from the beam axis. 

The upgraded UA2 calorimeter system |TB[ , shown in Fig. ^ covered the polar angular 
range, 6° < 9 < 174°, and was used to measure the central system, X. In order to isolate 
React. |I]from other (background) events, rapidity-gaps are imposed offline between p and 
X, and between p and X, by requiring the absence of charged-particle hits in the UA2 
Time-Of-Flight (TOF) counters. These counters are indicated in Fig. |6| and cover the 
range of pseudorapidity, 2.3 < |r/| < 4.1, in both arms (2°-12° and 168°-178°). Since the 
TOF counters have some overlap with the small-angle region of the end-cap calorimeters, 
the calorimeter minimum acceptance angle for the events considered here is increased from 
6° to 12° in both arms. 



2.1 Triggering 

Since the main goal of the UA8 experiment was to make measurements of hard- 
diffraction scattering in React. |2| [ffl H, UA8 was interfaced to the UA2 data acquisi- 
tion system, which allowed the formation of triggers based on various combinations of p 
and /or p momenta and transverse energy in the UA2 calorimeter system. Parallel triggers 
were also employed to yield samples of elastic and inelastic diffraction reactions with no 
conditions on the energy in the calorimeter system. 

In order to find evidence for React. [I], one of the supplementary triggers required 
detection of a non-collinear p and p pair. The p and p were both required to be either in 
the "UP" spectrometers (above the beam pipe), as shown in Fig. 0(a), or in the "DOWN" 
spectrometers (below the beam pipe). 

During the 1989 run, 1297 events were recorded in which both p and p tracks were 
detected and the calorimeter system had a total recorded energy greater than 0.25 GeV. 
The remainder of the event-selection procedure for these events is described in Sect. |3.1| . 

The essential topology characteristic of these events is summarized in Table |l|. It is 
seen that, when both p and p have x p > 0.95, 48% of the events have rapidity-gaps in 
both arms (with pseudorapidity, 2.3 < |?7| < 4.1). However, when one or the other of 
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the tracks has x p < 0.95, the percentage which possess rapidity-gaps in both arms falls 
to only a few percent. Thus, the first class of events appears to constitute a unique set, 
distinct from those events in which either p or p have x p < 0.95. 

A secondary data sample of React. was extracted from data which were triggered 
by requiring that either p or p is observed, as shown in Fig. 0(b). In these events, since 
there is no selection bias on momentum transfer, t, of the undetected particle, its natural 
distribution prevails, with an average value, \t\ ~ 0.2 GeV 2 . There were 62,627 such events 
recorded, after offline cuts for pile- up, beam-pipe geometry and halo cuts are made @. 
These data are dominated by single diffraction, React. |2|. However, we show in Sect. [T]2] 
that the offline imposition of the rapidity-gap veto in both arms isolates React. [1] in this 
data sample. 

3 Event selection 

We henceforth refer to the event sample for which both p and p were required in the 
trigger as the "AND" data sample. The events for which either p or p are detected are 
referred to as the "OR" data sample. 

3.1 "AND" data sample 

The four constraints of energy-momentum conservation can be examined in individual 
events because the entire final state of React. [I] is seen in the Roman-pot spectrometers 
and the calorimeter system. 

Fig. |8|(a) shows the distribution of total visible energy (p and p and central system, 
X) for these events. Although there is clearly a component of events which possess the 
full available energy of 630 GeV, a significant fraction of the events have less energy. 
Fig. ||(b) shows the same distribution, but for those events in which rapidity-gaps have 
been imposed using the Time-Of-Flight (TOF) counters in both arms. 188 events remain 
in a clean signal at 630 GeV. 

Having seen that the energy constraint is well satisfied, we now consider the three mo- 
mentum constraints. As implied by Fig. 0(a), a minimum accepted transverse momentum 
of ~ 1 GeV for each of p and p corresponds to a net transverse momentum imbalance, 
Pt > 2 GeV, which is compensated for by a corresponding (opposite) momentum vector 
in the UA2 calorimeter system. In order to observe this, we define a summed momentum 
vector in the calorimeter. The cell energies observed in the UA2 calorimeter system are 
summed up as (massless) vectors to approximate the total momentum vector, P(X), of 
the system X in React. [I]. The azimuthal angle of P{X), $x, is plotted in Fig. ^a) vs. 
the azimuthal angle of the summed momentum vector of the final-state p and p particles. 
There are peaks seen at 90° and 270°, corresponding to the cases where both p and p 
are both in their DOWN spectrometers or both in their UP spectrometers (as sketched 
in Fig. [0(b)), respectively. Although $x has no acceptance or trigger bias, in both cases 
it is seen to be opposite the azimuthal angle of the summed pp momentum vector in the 
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figure. 

The projection of the points in Fig. |S](a) on the $>x axis is shown in Fig. |3](b). The 
different intensities in the two peaks are due to small differences in the distances of the 
Roman-pots from the beam axis in the two spectrometers, resulting in a mismatch in 
their low— P t cutoffs. The solid curve is a Monte-Carlo calculation and shows that the 
width of the peaks is understood. We select 139 events with $x in the bands 90° ± 20° 
and 270° ± 20°. 

Although the summed transverse momentum of p andp, Pt(p+p), drops off sharply be- 
low 2 GeV, the transverse component of the calorimeter vector, P t (X), shown in Fig. |^(c) 
displays a much broader distribution due to the resolution of the calorimeter and the fact 
that, at small particle energies, some energy is lost before the particles reach the sensi- 
tive volume of the device. Fig. [10|(a) shows the transverse projection of the calorimeter 



momentum vector together with the result of a Monte-Carlo simulation p8| , P9f| of the 
UA2 calorimeter system. This shows that the UA2 calorimeter simulation software does 
a good job in describing the calorimeter's low energy response. We select 126 events with 
1 < Pt{X) < 3 GeV for further analysis. 

The degree of longitudinal momentum balance is demonstrated in Fig. |^(d), a his- 
togram of total longitudinal momentum, SP; „ 9 , which includes the p, p and calorimeter 
longitudinal energies. A final sample of 107 DPE events satisfies the selection, |£P; ong | < 7 
GeV. The shaded histogram in Fig. |8](b) shows the total visible energy for these events 
and demonstrates how well the energy constraint is satisfied. Fig. ^(c) is essentially the 
same as Fig. |3|(b), except that the order of the TOF veto and momentum-conservation 
cuts is inverted. 

Table H summarizes the event losses due to the cuts described here. They are compared 
with the effect of the same cuts on the Monte-Carlo generated events discussed in Sect. |. 
The similarity between the two sets of numbers implies that most of the 188 events shown 
in Figs. |8|(b) and ||(a,b) are in fact real examples of React [[]. 

An additional point can be made that there is an insignificant contribution in the data 
sample from events in which the observed proton comes from a diffractively-produced low- 
mass system (Baksay et al. measured that this occurs (12 ±2.5)% of the time). Such 
events would lead to an asymmetry and tail on the low side of the total visible energy 
distribution in Fig. || Although a small tail of this type does exist, it disappears when 
the momentum conservation cut is made. Thus we conclude that the rapidity-gap veto 
combined with momentum conservation eliminates this source of background. 



3.2 "OR" data sample 

As remarked above, the "OR" triggered sample is dominated by React. [| However, 
the small component which is React. [1] can be isolated by selecting those events which 
possess rapidity-gaps in both arms. 

A signature which distinguishes React. [I] from React. |2| is the presence of a longi- 
tudinally forward-backward symmetric distribution of particles in the UA2 calorimeter. 
Fig. |TT|(a) shows distribution of the summed longitudinal momentum component of all 
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struck cells in the calorimeter for the triggered "OR" data sample. In constructing this 
plot, each event is plotted on the negative side if the summed vector is in the same hemi- 
sphere as the observed trigger particle, and on the positive side if the vector is in the 
opposite hemisphere. A large asymmetry is seen favoring the hemisphere opposite the 
trigger particle, as expected for React. 0. 

Fig. |ll](b) shows the subset of events in Fig. |TT](a) which have no hits in any of the 
UA2 Time-Of-Flight counters. This corresponds to rapidity-gaps in the range 2.3 < 
I77I < 4.1 in both arms. The forward-backward asymmetry seen in the calorimeter system 
disappears. The "AND" events are also plotted in Fig. |TT|(b) and are seen to have the 
same summed calorimeter momentum distribution as do the "OR" events. We take the 
events in the resulting symmetric distribution as the candidates for React. [I]. 

The next step in the selection of React. [I] is to look at the equivalent of Fig. ||(b) 
for this sample, namely the azimuthal angle, $x, of the summed calorimeter vector. 
Fig. |12|(a) shows its distribution for those single-diffractive events in which the observed 
p or p is seen in the DOWN spectrometer. Not surprisingly, a correlation is seen between 
the azimuthal angles of the observed p or p and the summed calorimeter vector. However, 
when we make the "OR" event selection, by imposing the rapidity-gap condition using 
the TOF counters, we see in Fig. 0(b) that the correlation becomes much stronger. The 
distribution is broader than that seen in Fig. |9](b) for the "AND" events because of the 
unknown (small) P t of the unobserved final-state p or p and also because of the smaller 
energy in the calorimeter. The Monte-Carlo simulation of the UA2 calorimeter is in 
reasonable agreement with the observed distribution. We select 698 "OR" events, with 
$x either in the range 90° ± 20° or 270° ± 20°. 

Finally, for these "OR" candidates, we examine the summed transverse momentum in 
the calorimeter shown in Fig. |10](b) . Because this vector is opposite only one observed 
vector of the p or p, its average value is less than that seen in Fig. [13(a) for the "AND" 
sample. However, it also is in reasonable agreement with the Monte-Carlo simulation of 
the UA2 calorimeter. 

Figure [TB| shows that the momentum transfer (t) distribution of the "OR" events is 
in good agreement with that of the full single-diffractive data sample. This is consistent 
with our basic assumption that the flux factor is common to Reacts. |2| and [|. The 
lower statistics "AND" data sample (not shown here) is also compatible with the single- 
diffractive data. 



3.3 Feynman-Xp distribution 

The shaded distributions in Figs. [L4|(a,b) show the distributions of Feynman-x p and 
Xp for the final "AND" and "OR" data samples, respectively. They are essentially indis- 
tinguishable. The open histogram superimposed on both "AND" and "OR" distributions 
(shaded) in Fig. [TJ] is the x p /xp distribution in the single-diffractive data of React. ^ 
in our experiment 0. In order that both sets of distributions have the same kinematic 
conditions, the single-diffractive data are plotted only for those events that have no hits 
in the TOF counters on the trigger side, which cover pseudorapidity, 2.3 < 77 < 4.1. Each 
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open histogram is normalized to its shaded distribution for the bin: 0.990 < x p < 0.995. 

We see that the single-diffractive data possess a significant event population for x p > 
0.995, which is not seen in either set of data for React. [l|. As discussed in the Introduction 
and in Sect. ^, this apparent breakdown of factorization is merely a kinematic suppression 
in React. |1|, due to the requirement that s' > for the two-Pomeron system. In the 
Monte-Carlo generation of two independently-emitted Pomerons according to Eq. [| (see 
Sect. ||), 61% of all events, in which both p and p have \t\ > 1 GeV 2 , have s' < and are 
discarded. The rejected events mostly have small-£; for example, when either Pomeron 
has £ < 0.0005, 100% of the events are rejected, whereas when £ « 0.03, only 26% of 
the events are rejected. This qualitatively accounts for the difference between DPE and 
single diffractive data near x p = 1 in Fig. [TJ|. For the "OR" topology, when only one 



\t\ > 1 GeV 2 , the rejection at small-£ is about the same (95%), whereas there is more 
rejection at £ = 0.03 (51%); this accounts for the small difference in shape between "AND" 
and "OR" data in Fig. [TJ|. The detailed shapes of these distributions depend on cr^-p, 
which we have not yet determined. 



4 Monte— Carlo event generation 

A complete Monte-Carlo simulation [^] of React. [TJ was performed to determine the 
spectrometer and calorimeter acceptances as well as the efficiencies of the various cuts. 
Events were generated such that the Pomerons are emitted independently from proton 
and antiproton, respectively, according to Eq. |5|, using the Pomeron flux factor [|1^, [J 
in Eq. [| afej, is assumed to be independent of s', although in Sect. || we will look for 
departures from this assumption. 

Points were chosen randomly in the 6-dimensional space[], (t 1 , £ 1; 4> 1 , t 2 , £2, 4>2), ac- 
cording to the product of two flux factors. Each such point defines the properties of the 
Pomeron-Pomeron system, its energy and its momentum vector. We have observed that, 
even though the two Pomerons are assumed to be independently emitted, not all points 
in the 6-dimensional space are kinematically allowed because the associated Pomeron- 
Pomeron invariant mass may be unphysical (i.e., s' < 0). Thus, events are retained only 
if they are in regions of the 6-dimensional space for which s' > 0. In our |t|-domain, 
1-2 GeV 2 , such events are 39% of the total generated. We note that, even though the 
Pomerons are generated isotropically and independently in azimuthal angle, 0, correla- 
tions occur due to this kinematic suppression. 

The number of particles of the central system, X, in React. [I] is generated according 
to a Poisson distribution with its mean charged particle multiplicity depending on Mx, as 
measured in a study of low-mass diffractive systems [[JTJ, n = O.QMx (Mx = Vs* in GeV); 



the mean number of neutral particles is assumed to be one-half the number of charged 
particles. The tracks are generated isotropically in the Mx center-of-mass (see Sect. |5|). 



As described in Sect. |572|, where the central system is seen to have longitudinal structure 



5 Since the 3 observables of a final-state proton or antiproton are uniquely related to those of its 
associated "Pomcron, we use the 'Pomeron variables, £, t and 4>. 
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for Mx > 5 GeV, the Monte-Carlo generator is tuned to agree with the data. 

After phase-space generation of the complete events, their data were passed through 
detector simulation software for both the UA8 spectrometers and the UA2 detectors [28] , 
and then through the same offline analysis software and cuts used for the real data. 

As already discussed in Sect. [| Table ^ shows the good agreement between the real 
event losses with the momentum cuts described in Sect. |]and those on the Monte-Carlo 
events calculated here. Since the Monte-Carlo event sample of React. [I] suffers a 46% loss 
when the TOF veto is imposed, the net efficiency for event retention due to TOF veto 
and event selection cuts is 26%. 

The combined geometric and detection efficiency of proton and antiproton is about 
6 • 10~ 4 at an average \t\ of 1.2 GeV 2 . Figures [T3|(a,b) show the M^-dependence of the 
overall geometric and detection efficiencies averaged over all other variables, for "AND" 
and "OR" data, respectively, when the observed particles are in the range, 1.0 < —t < 
2.0 GeV 2 . The 26% central system detection efficiency is also included in these efficiencies. 
The fall-off in acceptance for the "AND" data at low mass results from the fact that, 
kinematically, low-mass events tend to have back-to-back p and p, which do not satisfy 
the "AND" trigger topology seen in Fig. [<]. The "OR" trigger topology does not have 
such a bias against low-mass events. 



The longitudinal structure seen in Sect. |5.2j for central system masses larger than 
about 5 GeV must impact the acceptance of the central system. Thus, we show the 
"AND" acceptances in Fig. |15|(a), first assuming isotropic decay of the central system 
and then, for Mx > 4 GeV, with a longitudinal decay distribution which matches that 
which is observed. At large mass, the acceptance for longitudinal decay is 25% smaller 
than the acceptance for isotropic decay. Since, with our statistics, it is not possible to 
properly study the transition from isotropic to longitudinal decay, we take the effect into 
account in the following way. For Mx < 4 GeV, we use the calculated acceptance for 
isotropic decay, shown as the solid curve. For Mx > 10 GeV, we use the fitted horizontal 
solid line in the figure. In the intermediate range, between 4 and 10 GeV, the dotted 
interpolation line is used. 

For the "OR" acceptance at low mass (M x < 4 GeV), we use the calculated acceptance 
for isotropic decay in Fig. [l5|(b), just as we did for the "AND" data. For Mx > 10 GeV, 
we use the solid horizontal line which is 25% below the calculated efficiency for isotropic 
decay. Again, the dashed interpolation line is used in the intermediate region. 



5 Calorimeter measurement of central system 

We use the UA2 calorimeter information to study the invariant mass and other prop- 
erties of the central system, X, in React. |l[ The UA2 detector simulation software [p£ 
was used to perform a complete Monte-Carlo study of the calorimeter response. As noted 
above, the UA2 simulation software is remarkably good in describing the low-energy de- 
posits encountered in our data. 



10 



5.1 Invariant mass of central system 

Since we do not directly observe the individual particles of this system, but rather the 
energies deposited in the calorimeter cells, we assume that the non-zero energy in each 
"struck" cell of the calorimeter is caused by a massless particle, and then calculate, 

Mobserved 2 = (^Ei) 2 — |EPj| 2 , (7) 

summing over all cells. 

Figure [16] is the result of a Monte-Carlo study |29| which shows that this procedure 



underestimates the true mass by an amount that increases with mass. This effect results 
from incomplete detection of energy; for example, the finite cell size leads to overlapping 
energy deposits from neighboring particles, or some energy can be lost before the particles 
enter the calorimeter. The difference between the true MC mass, M trU e, and the calculated 
or observed mass, M b ser ved, is plotted vs. M true . The dependence is well fit by the 
equation (with M in units of GeV): M true — M b se rved = (1 + 4M true )/14, which can be 
rewritten as: 

M true = lAM observed + 0.1. (8) 

We define the corrected mass, M\ = s', to be the true mass given by this equation and 
only refer to these corrected values in the remainder of this paper. 

The validity of the calorimeter invariant mass calculation may be conveniently tested 
by comparing it with the missing mass calculated using the measured p and p 4- vectors for 
an event. Although the experimental uncertainty in a "missing mass" calculation is much 



larger than for the calorimeter invariant mass, they should agree on average. Fig. [T7| shows 
the average missing mass calculated for the events in each of the calorimeter invariant mass 
bins shown in the figure. The observed clustering of the points around the diagonal and 
the absence of any systematic shifts constitutes proof that the calorimeter mass evaluation 
is reliable. 

The Mx distributions of the system X in React. [I] are shown for the final selected 
"AND" and "OR" event samples in Figs. |T5|(a,b), after requiring that the momentum 
transfer of all detected protons and antiprotons be in the range, 1-2 GeV 2 . From the 
relatively flat acceptance curves in Figs. |T3|(a,b), we see that the observed shapes of the 
distributions are reasonably good representations of the true distributions (except for the 
lowest mass bin in the "AND" data) . In Sect. FO, we show that part of the low-mass 



enhancements are attributable to an explicit s' dependence in crtpp, corresponding to an 
enhanced Pomeron-Pomeron interaction in the few-GeV mass region. However, with an 
estimated 1.8 GeV mass resolution obtainable from the calorimeter, we are unable to 
observe details of any possible s-channel resonant structure in this spectrum. 

5.2 Other properties of the central system 

In addition to the invariant mass distribution of the system, X, other properties of the 



system can be studied. One is the particle multiplicity of the central system. Figure |19 



shows the number of calorimeter cells struck as a function of the corrected calorimeter 
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mass, Mx- The solid line is a fit to the data; the dashed line is based on the naive 
multiplicity expectation assuming < iV >= 0.6Mx (Mx in GeV) for the number of 
charged particles (7r + and tt~). The number of 7r° is assumed to be Poisson-distributed 
with a mean of 0.3M X ; each of these is assumed to appear as two 7. The resulting dashed 
line is the function iV = 1.2Mx and clearly captures the gross features of the data. The 
observed numbers of struck cells lie somewhat above the line, as expected geometrically 
from the cluster widths in the calorimeter and the finite cell sizes. A complete Monte- 



Carlo simulation p{| accounts for the small observed differences. We can conclude that 
the number of observed struck cells increases with mass roughly as expected, and the 
total observed multiplicity displays no anomalous features. 

Because of the separated electromagnetic and hadronic section of the UA2 calorimeter, 
it has also been possible to study the fraction of electromagnetic energy possessed by the 
central system. Fig. ^0] shows the distribution in this fraction for the "AND" events in 
the low-mass enhancement, Mx < 6 GeV, compared with the Monte-Carlo generated 
distribution assuming, on the average, equal numbers of tt + , tt~ and tt° in the track 
generation. Again we see no anomalous features in this variable. The enhancement visible 
in both data and Monte-Carlo when the ratio, (e.m. energy) /(total energy), equals unity 
corresponds to low mass systems where the slow pions deposit all their energy in the 
electromagnetic calorimeter cells. 

We have examined the angular distributions of calorimeter cell energies in the center- 
of-mass of the X system, dN/dcosd, with respect to the Pomeron-Pomeron direction of 
motion. Figs. |2T](a,b) show these distributions for Mx < 5 GeV and Mx > 5 GeV. At 
the higher masses we see a similar type of forward-backward peaking as is seen in all 
hadronic interactions as a result of the presence of spectator partons. In the present case, 
this would imply that there are spectator partons in the Pomeron. We have already 
reported 0, [3^] similar effects in Pomeron-proton interactions in the single- diffract ive, 
React. 0. 

The Monte-Carlo histogram in Fig. |2l|(a) shows isotropically-decaying events. In 
Fig. |2T|(b), the histogram shows a Monte Carlo event sample which has been selected in 
such a way that it has the same forward-backward peaking as the experimental distri- 
bution. For each isotropically-decaying Monte-Carlo event, the mean value of cos 2 9 is 
evaluated averaging over all outgoing tracks in the central system. We have found 
that, if Monte-Carlo events are selected for which this quantity is larger than 0.375, the 
selected events follow the experimental cosO distribution. 



6 Cross sections 

The observed mass distributions, dN/dMx, shown in Figs. |18|(a,b) for the "AND" 
and "OR" data, respectively, are converted into absolute cross section distributions, 
daupE/dMxi in the following way. Bin-by-bin, the numbers of events in Figs. are 
divided by the Monte-Carlo acceptance curves in Figs. [I5|(a,b). Then, all are divided by 



a global efficiency, eo, for event retention when halo and pileup cuts [12, 29] are made, 
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and by the appropriate integrated luminosity for each trigger sample: 



The "AND" ("OR") data samples have an efficiency, eo, of 0.54 (0.76) and an inte- 
grated luminosity / Cdt, of 2894 ntr 1 (5.4 ntr 1 ). The "OR" luminosity is "effective", 
due to prescaling of the "OR" trigger. In both cases, the cross section is given only for 
momentum transfer of the observed trigger particle(s) in the range: 1.0 < \t\ < 2.0 GeV 2 . 
In the "OR" case, the unseen p or p has its "natural" P t distribution and therefore peaks 
at small values. Thus, the observed cross section is much larger for the "OR" data. The 
resulting cross sections for the two triggered data samples of Reaction [T], da/dM x , are 
the points shown in Figs. |22"|(a,b). 



6.1 'Pomeron— Pomeron total cross section 

We now extract the Pomeron-Pomeron total cross section, <7pp, from the data, so 
that we can look for deviations from our earlier assumption that it is independent of Mx- 
The histograms in Figs. ^(a,b) are Monte-Carlo predictions for the daopE/dMx points 
in the figure, made using Eq. ^|and Fp/ P (t,£) in Eq. [| We assume a constant Pomeron- 
Pomeron total cross section of crpp = 1 mb and the (arbitrary) value, K = 0.74 GeV -2 , 
as discussed in Sect. [L| 

Since the UA8 spectrometer acceptances || do not vary significantly over the range of 
£ studied here, the ratios of the points to the histogram values in Figs. |22| give values of the 
Pomeron-Pomeron total cross section, Opp, vs. Mx- These ratios are shown in Fig. [23| 
for both the "AND" and "OR" data. We note that, despite the large difference between 
the measured cross sections, daopE/dMx, in Figs. |2|(a,b), both data sets yield the same 
general properties for Opp: enhancements for Mx < 8 GeV and relatively Mx independent 
shapes at larger Mx- Although the two values of crpp in the first bin (Mx < 2 GeV) are 
consistent with being equal, in the next three bins (2 < Mx < 8 GeV) the "AND" cross 
sections are about three times larger than the "OR" cross sections. However, only the 
statistical errors on Opp are shown in Fig. ^| We discuss their systematic uncertainties 



in the following section. 

The small-Mx enhancement in Fig. ^ also reflects itself in the observed x p (x p ) 



distributions seen above in Figs. TA. Such a correlation must exist because of the kinematic 
relation, Mx 2 = s' = £i^2s; small Mx correlates with small £. Fig. [24] repeats the x p (x. 



pj 



distribution for the "AND" data in Fig. |T4](a). The solid curve normalized to its area is the 
Monte-Carlo prediction which assumes an s'-independent ffp-p. The pronounced excess 
of events near x p = 1.0 in the experimental distribution, compared with the Monte-Carlo 
distribution, is another manifestation of the low-mass enhancement in Opp. 

The low-mass enhancements seen in both distributions in Figs. |2^ are most likely too 



large [33] to be due to a breakdown of factorization at small mass, especially for the "AND" 
data. Thus, the rise may indicate that glueball production is a significant component of 
the low-mass Pomeron-Pomeron interaction, although not necessarily an s-channel effect. 
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That is, the observed invariant mass could be that of a glueball plus other particles, which 
would not lead to resonance structure in the mass distribution. In any case, with a mass 
resolution of ~ 1.8 GeV, no s-channel structure could be seen. 



6.2 Systematic uncertainties 

We now ask how possible systematic uncertainties influence the interpretation of the 
a-pp results shown in Fig. A systematic uncertainty common to both sets of cross 
section results comes from the Monte-Carlo acceptance shown in Fig. |Tj| We estimate 
that this uncertainty is smaller than 15%. In the case of the "AND" results, the statistical 
errors are much larger than this. 

There is an additional systematic uncertainty that is very different for the "AND" and 
"OR" results. This arises from the possible non-universality of the Pomeron flux factor. 
It is already known that F-p/ p (t, £) is not universal between ep and pp collisions because of 
the different effective Pomeron trajectory intercepts found in the two cases, attributable 
to different multi-Pomeron-exchange effects. Similarly, if multi-Pomeron-exchange is not 
identical in Reacts. [3] and |2|, there would be some uncertainty as to whether the same flux 
factor should appear in both Eqs. ^] and |3|. 

We note, however, that this potential uncertainty does not exist for the "AND" results, 
because both final-state baryons have \t\ > 1.0 GeV 2 , where all the evidence flT7|, [HJ points 
to an s-independent Pomeron trajectory; in that high-|t| region, F-p/ p (t,£) appears to 
be insensitive to the damping which mainly leads to an s-dependent effective Pomeron 



intercept [Ty] at t = 0. Thus, the statistical errors on the "AND" results in Fig. ^3 
dominate any uncertainties from this source. 

The situation is very different for the "OR" data sample because the unseen final- 
state baryon has low-|i| and its flux factor in Eq. [| is sensitive to the choice of effective 
e value. We have thus recalculated the "OR" cross sections in Fig. ^ assuming a larger 
Pomeron trajectory intercept, e = 0.10, in F-p/ p (t,£). This is an extreme (unrealistic) case 
which assumes that there are no damping contributions from multi-Pomeron-exchange in 
React. |] at our energy. We find that the "OR" cross sections decrease from those shown 
in the figure; for example, the lowest mass point decreases by 58%, whereas the point at 
11 GeV decreases by 41%. Thus, any increase in the effective e used in calculating cr^-p 



from the "OR" data increases the disagreement already observed in Fig. ^3. 

The final systematic uncertainty comes from our model assumption that there are 
no azimuthal angle correlations between final-state p and p in Eq. |5|, other than the 
kinematic one referred to in Sect. [I]. However, the following two facts lead us to suspect 
that this may not be true: (a) there are differences between the "AND" and "OR" 
a-pp results, and (b) the "AND" and "OR" data samples have very different Pomeron- 
Pomeron configurations. The difference in Pomeron-Pomeron configurations is most 
easily visualized with reference to our Fig. [!]. In the upper figure ("AND"), it is seen that 
both Pomeron transverse momenta, P t , are approximately in the same direction; hence 
AP t fa 0. In the lower figure ("OR"), one transverse momentum is near zero and the 
other is near 1 GeV; hence the "OR" data sample corresponds to AP t ~ 1.0 GeV. 
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Thus, the observed differences in a% for the "AND" and "OR" data in Fig. || at low 
mass suggest that a^p depends on the Pomeron-Pomeron relative configuration and is 
larger at low mass when the two Pomerons move approximately in the same direction. 
We also note that the systematic uncertainty due to our limited understanding of multi- 
Pomeron-exchange effects can not be responsible for the effect now being discused, since 
it can only make the difference larger. We return to the physics of the present discussion 
in the concluding Sect. |7|. 



6.3 Test of factorization 



We now test the factorization relation, Eq. |6], between a^-p, aip l p and a 
both sides of Eq. | by K 2 , we find: 



tot 

PP ' 



K 2 



a 



tot I 



<T*; iaZ (s') 



If we multiply 



(10) 



which relates precisely the measured quantities. Thus, it is evident that tests of factor- 
ization using Eq. |6| or [T^ are equivalent and we may, with no loss of precision, assume the 
value, K = 0.74 GeV -2 (see discussion in Sect. P, and make the test using Eq. §. 

In order to calculate the right-hand-side of Eq. || a function of s', we use the following 
parametrizations for the Pomeron-proton total cross section ||: 

°' 72 -[( S r io +4.o( S ') 



a 



tot 
Vp 



-0.321 



0.74 



mb, 
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and for the proton-proton total cross section |3J, [35] : 

27 • <r a50 + 55 



a 



tot 

PP 



,0.10 



-0.32 



mb. 



(12) 



These functions are shown as the dashed curves in Figs. |25|(a,b). Since Eq. || is only valid 
for the Pomeron-exchange component of these functions, we show only the first terms in 



Eqs. O and 12 as the solid curves in the figures. 



The dashed line in Fig. |23] shows the factorization prediction for cr^p calculated using 



Eq. p and the Pomeron terms in Eqs. 11 and 



a. 



tot 

VP 



V) 



(j 



tof\2 
Vp\ 



a 



total 
PP 



(0.72/0.74) 2 
18 



J\0.10 



(13) 



We see that there is increasingly better agreement between the prediction and the mea- 
sured crpp points as the mass increases, as expected, since the measured points contain 
both Pomeron exchange and P.eggeon exchange. The results are seen to be in reasonable 
agreement with the validity of factorization for Pomeron-exchange in these reactions. 

The solid curve in Fig. |23| is a fit to the "OR" points with Mx > 10 GeV of the sum 
of Eq. [1^ and a P.eggeon-exchange component as in Eq. H : 

(0.72/0. 74) 2 



/\0.10 



+ R-(s' 



-0.32i 



(14) 



We find a value, R = 13.6 ± 4.7, with a x 2 /D.F. = 1.3. 



15 



7 Conclusions and predictions 



Table demonstrates that there is a new class of events, the so-called double-Pomeron- 
exchange (DPE) events with characteristic rapidity-gaps in both arms, which appear when 
both x p values are greater than 0.95. The analysis which follows this observation shows 
that, remarkably, the Regge formalism describes all inclusive DPE and single-diffractive 



data, using the empirical s-dependent effective Pomeron trajectory of Ref. [fL9|| , which 
is believed to be due to increasing multi-Pomeron-exchange effects with energy |TH . 
That Regge phenomenology works as well as it does, despite the complications of multi- 
Pomeron-exchange, should place constraints on a theory of such multi-Pomeron-exchange 
effects yet to be developed. 

We observe that the produced central systems in the DPE events display no anoma- 
lous multiplicity distributions or electromagnetic energy fraction of the total observed 
calorimeter energy. Our measurements agree with normal expectations. 

The main result of the work reported here is the Pomeron-Pomeron total cross sec- 
tion, cjpp, shown in Fig. At large mass, there is a statistically weak agreement with 
factorization predictions. However, for Mx < 8 GeV, both the "AND" and "OR" results 
exhibit large enhancements in u^p, with the "AND" result being about three-times larger 
than the "OR" result. There may be a dynamical reason for this, based on the fact that 
the "AND" and "OR" data sets have different Pomeron-Pomeron configurations. 

These results are probably related to the WA102 Collaboration observations |^| , at 
the much lower center-of-mass energy^, y/s = 29 GeV, that the production of known 
quark- ant iquark states and glueball candidates depend differently on the difference in 
transverse momenta of the two Pomerons in React. [I]. The primary WA102 effect is that 
the production of qq states appears to vanish as AP t — > 0. Close and Kirk |37| suggested 
that this effect may serve as a "glueball filter" . On the basis of the WA102 results, Close 
and Schuler [38] argue that the effective spin of the Pomeron can not be zero and that 
the Pomeron transforms as a non-conserved vector current. 

Taken at face value, the WA102 results [^] imply that our "AND" data do not contain 
any qq states and, hence, that the enhancement for 2 < Mx < 8 GeV in those data may 
be due to production of some glueball-like objects. There would be a mix of s-channel 
production (i.e. glueball alone) and production with other particles. The observation of 
resonant mass structure is precluded in the present experiment because of our poor mass 
resolution of ~ 1.8 GeV. 

The next generation experiment of the type reported here should utilize a central 
detector capable of detailed studies of the produced central systems (including particle 
identification). In order to be able to observe the azimuthal angle correlations of the 
outgoing protons with properties of the central systems, the Roman-pot systems (on 
both arms) should have as full an azimuthal coverage as possible. 



6 This translates to "Pomeron momentum fractions, £, which are 22 times larger in WA102 than in 
UA8. Since this can lead to non-Pomeron-exchange contributions [§, as large as 50% in WA102, we 
should not expect to observe exactly the same effects. 
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7.1 Predictions for Tevatron and LHC colliders 



The fact that the "AND" and "OR" data yield essentially the same er-pp at larger mass 
implies that our rapidity-gap procedures for isolating React. |l] are fundamentally sound 
and can used at the much higher energy experiments at the Tevatron and LHC to look 
for rare states with greatest sensitivity. 

The study of the relatively pure gluonic collisions in React. [3] at higher energy colliders 
may yield surprising new physics. In the UA8 papers on the occurrence and study of jet 
events in single-diffraction 0, it was reported that, in about 30% of the 2-jet events, 
the Pomeron appeared to interact as a single hard gluon with the full momentum of the 
Pomeron (the so-called "Super-Hard" Pomeron) . This result suggests that, in roughly 
10% of hard Pomeron-Pomeron interactions, there is effectively a gluon-gluon collision 
with the full Mx of the central system. Thus, for example, at the LHC with \fs = 14 TeV, 
a rather pure sample of central gluon-gluon collisions should occur with Mx as large as 
0.03(14) = 0.42 TeV (remember from Sect. [I] that we believe there is essentially pure 
Pomeron exchange at £ = 0.03). 

The phenomenology developed in our study of React. [I] can be used to make cross sec- 



tion predictions at the Tevatron (y/s = 2 TeV) and at the LHC (-y/i = 14 TeV). Figure |26 



and Table |3] show the results of Monte-Carlo calculations of da^PE/ 'dM x (integrated 
over all t), assuming that crpp is M^-independent and constant at 1 mb. Since the fitted 
value |19| of the effective Pomeron-Regge-trajectory intercept, 1 + e = 1.035, at y/s = 
630 GeV, is also compatible with the available data at the Tevatron, we give the results 
for e = 0.035 at both Tevatron and LHC. Schuler and Sjostrand suggest, in a model 



of hadronic diffractive cross sections at the highest energies, that e = is a reasonable 
approximation and we therefore also give results for this value. The observed peaking at 
small mass directly reflects the ^-dependence of the Pomeron flux factor in the proton. 

To obtain cross section predictions for central Higgs production ||0| in React. [I] from 
the "Super-Hard" component of the Pomerons , the calculations in Fig. |26] can be mul- 
tiplied by the calculated QCD cross section for the process (in units of mb) based on the 
best available Pomeron structure function. 

We note that the cross section predictions obtained in this way give the total cross 
section for React. [I], where there are no selection cuts on either final-state p or p. Clearly 
this yields the largest sensitivity for rare events. In that case, if rapidity-gap vetos are 
used to suppress background events, corrections must be made for the acceptance loss of 
central system particles in the rapidity-gap regions. 



7.2 Predictions for forward spectrometers 

For completeness, it may be useful to briefly summarize the possibilities for detection 
of DPE processes using existing or planned forward multiparticle spectrometers. There 
are two classes of such experiments, those traditionally called fixed-target experiments 
and those installed at storage- ring colliders [CT. 

Forward spectrometers installed at colliders can observe DPE processes if there is 
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an asymmetry between £j and £2- The earliest example of such a measurement was 
Experiment R608 studying pp interactions with = 63 GeV at the Cern Intersecting- 
Storage-Rings. Central production of .0(1285) was observed ||42|| with almost pure helicity 
±1, later explained by Close and Schuler as due to the Pomeron behaving as a non- 
conserved vector current. In that process, the Pomeron appeared to dominate even though 
£i — £2 ~ 0.35. In the future, the higher energy forward-spectrometer B-experiments, 
LHCb and B-TeV, will have good access to low-mass central systems with much smaller 



values of £, as was pointed out in the LHCb Letter-of-Intent |43f . 

In fixed-target experiments, centrally-produced systems are boosted forward with the 
7 of the center-of-mass, such that Eiab = jMx- Experiment WA102 |36| was the first 
experiment to carry out major DPE studies using this approach although, as commented 
above, with a beam energy of 450 GeV, the Pomeron ^-values were larger than desired. 
We note that the existing experiment, HERA-BjyJ, running at the HERA 920 GeV 
proton storage ring using wire targets, could improve on the WA102 measurements. For 
example, with y/s = 42 GeV, production of a 2-GeV central system occurs with an average 
£ = 0.047. Fig. ^ shows daopEl 'dMx for pp interactions with a beam energy of 920 GeV, 
calculated as for Fig. |26|. As indicated after Eq. |], the Pomeron trajectory intercept [pTOf] 
used at this cm. energy is: 1 + e = 1.087. 
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p or p has x p > 0.95 
Other has x p in bin 


Number 
Events 


Fraction with 
rapidity-gaps 


0.70-0.75 


42 


0.02 ±0.02 


0.75-0.80 


113 


0.04 ±0.02 


0.80-0.85 


147 


0.07 ±0.02 


0.85-0.90 


163 


0.03 ±0.02 


0.90-0.95 


219 


0.06 ±0.02 


0.95-1.00 


314 


0.48 ± 0.04 



Table 1: Numbers of events and their fractions which have rapidity-gaps, 2.3 < r] < 4.1 
in both arms, for different x p selections of p and p. 1297 events have two reconstructed 
tracks and at least 250 MeV of energy in the calorimeter system. The table shows the 
998 events, in which either p or p has x p > 0.95, while the other has x p in the indicated 
bin. In the remaining 299 events, for which both p and p have 0.70 < x p < 0.95, only 
0.7% possess both rapidity-gaps. 



After 


Events 


percentage 


MC percentage 


Cut 


remaining 


remaining 


remaining 


TOF veto 


188 






$ x 


139 


74 ± 6 % 


72 % 


Pt 


126 


67 ± 6 % 


60 % 


y p 

z -' 1 long 


107 


57 ± 5 % 


50 % 



Table 2: Comparison of event losses for data and Monte-Carlo events (after TOF veto) 
as a function of the 3-momentum cuts in the event selection. 
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M x 


da DPE /dM x (mb/GeV) 




Tevatron 


LHC 


(GeV) 


e = 0.00 


e = 0.035 


e = 0.00 


e = 0.035 


1 


1.25E-01 


2.69E-01 


2.07E-01 


7.55E-01 


3 


5.44E-02 


1.00E-01 


9.42E-02 


2.96E-01 


5 


3.01E-02 


5.17E-02 


5.39E-02 


1.58E-01 


10 


1.31E-02 


2.11E-02 


2.53E-02 


6.76E-02 


20 


5.33E-03 


7.65E-03 


1.15E-02 


2.75E-02 


30 


3.02E-03 


4.14E-03 


7.17E-03 


1.62E-02 


40 


1.96E-03 


2.57E-03 


5.08E-03 


1.13E-02 


50 


1.40E-03 


1.78E-03 


3.86E-03 


8.34E-03 


60 


1.01E-03 


1.25E-03 


3.11E-03 


6.45E-03 


70 


7.76E-04 


9.42E-04 


2.56E-03 


5.23E-03 


80 


6.08E-04 


7.15E-04 


2.17E-03 


4.32E-03 


90 


4.78E-04 


5.54E-04 


1.85E-03 


3.75E-03 


100 


3.72E-04 


4.21E-04 


1.64E-03 


3.20E-03 


120 


2.36E-04 


2.48E-04 


1.28E-03 


2.46E-03 


140 


1.37E-04 


1.61E-04 


1.04E-03 


2.01E-03 


160 


7.45E-05 


9.35E-05 


8.87E-04 


1.60E-03 


180 


3.10E-05 


3.74E-05 


7.58E-04 


1.36E-03 


200 


1.12E-05 


4.21E-06 


6.55E-04 


1.18E-03 



Table 3: Predictions for daopE/dMx (mb/GeV) at Tevatron and LHC assuming an 
Mx-independent a^j, = 1 mb, for two values of effective Pomeron-trajectory intercept 
(see Fig. |(| 
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Figure 1: (a) Inclusive double- Pomeron-exchange reaction (the central blob is the 
Pomeron-Pomeron interaction) and its corresponding Triple- Regge diagram in (b); (c,d) 
Inclusive single diffractive reaction and its corresponding Triple-Regge diagram. In both 
cases, Pomeron-exchange dominance means % — j — "Pomeron. k can be either Pomeron 
or 7?.eggeon. 
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Figure 2: Raw (uncorrected) Feynman-x p distributions for different bins of momentum 
transfer (units are GeV 2 ) in single-diffractive data from the UA8 experiment JJ|. As 
explained in the text, the x p -dependence of the geometrical acceptance is not responsible 
for the observed peaks. 
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Figure 3: Ratios of forward elastic amplitudes which are equal if factorization is valid. 
The optical theorem implies that the following relation between the corresponding total 
cross sections should be valid if Pomeron-exchange dominates: cr^p = (^p) 2 / "™- 
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Figure 4: Particle trajectories through a UA8 Roman-pot spectrometer. The labels, 
"Quads" , refer to the low-/? machine quadrupole magnets. The center of the UA2 detector 
is at z = at the left of the sketch. The vertical lines indicate the positions of the UA8 wire 
chambers in the Roman-pots. The solid curves show the trajectories of 300 GeV particles 
(x p ~ 0.95), as described in the text. The shaded area shows the allowed trajectories for 
x v = 1 tracks. 
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Figure 5: UA8 spectrometer aperture viewed from the interaction region. The shaded 
rectangles indicate the sensitive regions of the first wire chambers at a distance z = 13 m 
from the interaction region center. The curved line indicates the walls of the beam pipe 
inside the quadrupole magnets. 
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CALORIMETERS 




Figure 6: A cross sectional view of the upgraded UA2 apparatus. Detectors which 
were used for the measurements reported here are the Calorimeters, the Time-Of-Flight 
(TOF) counters and the Silicon Vertex Detector within the Central Detector assembly. 
The TOF counters covered pseudorapidity from 2.3 to 4.1 in each arm. 
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Figure 7: (a) Side view sketch of an accepted event in which both p and p go into the UP 
spectrometers. In this case, the central system recoils downward with a minimum P t of 
~ 2 GeV. (b) Sketch of an event triggered on p or p. In this case, the unobserved p (or p) 
has momentum transfer close to zero on average and the central system recoils downward 
with a minimum P t of ~ 1 GeV. 
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Figure 8: (a) Total visible energy in "AND" triggered events (1297) with total calorime- 
ter energy, HE > 250 MeV, selected offline; (b) Open histogram contains events (188) 
after veto using TOF counters, as discussed in the text. Shaded events (107) are after 
momentum conservation cuts (Figs. |9|); (c) Open histogram are events (193) after mo- 
mentum conservation cuts, but before TOF cuts. Shaded events are after the TOF cuts 
(107 events). 
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Figure 9: "AND" events (188) after TOF veto selection (see text); (a) Scatter plot 
of the azimuthal angles of the p + p system and the summed calorimeter momentum 
vector; (b) $x projection of (a). The curve on the right-hand peak is a Monte-Carlo 
simulation described in the text; (c) \Pt(p + p)\ vs. Pt(X) measured in calorimeter, for 
events (139) which satisfy $ x selection in the bands, 90° ± 20° and 270° ± 20°. (d) 
Summed longitudinal momentum, SP; oriS , of p, p, and calorimeter for events (126) which 

— * 

satisfy the selection, 1 < P t (X) < 3 GeV. 107 events satisfy the cut, \Y>Pi ong \ < 7 GeV. 
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Figure 10: Transverse momentum, P t , measured in the UA2 calorimeter system: (a) 
"AND" data sample, projection of Fig. ||(c) (139 events ); (b) "OR" data sample (698 
events, after TOF and $x cuts). The histograms are from a Monte-Carlo simulation 
of the UA2 calorimeter response. 
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Figure 11: For "OR" triggered events, comparison of the summed longitudinal momen- 
tum in the calorimeter, with and without rapidity-gap veto (TOF). (a) Without TOF 
veto (the single-diffractive data sample); events are plotted on the positive axis if their 
summed momentum is in the hemisphere opposite the observed trigger particle. Events 
are required to have at least 250 MeV energy in the UA2 calorimeter system (partial 
sample, 15,080 events); (b) With rapidity-gap veto (1985 events). Solid points show the 
107 "AND" events, normalized to the "OR" data. 
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Figure 12: Calorimeter azimuthal angle, $x- (a) single-diffractive data sample, with p 
or p in DOWN spectrometer (&p + p selection in the band, 270° ± 20°) and with total 
calorimeter energy, HE > 250 MeV (5547 events); (b) Same as (a), but after TOF veto 
selection to obtain (partial) "OR" data sample (635 events). The curve in (b) is the 
Monte-Carlo simulation described in the text. 
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Figure 13: Observed momentum-transfer distributions for the "OR" data sample of 
React. [I] (solid points). The histogram normalized to the points is the single-diffractive 
data, React. |[ 
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Figure 14: Feynman-x p distribution; (a) Shaded histogram is x p and x p (2 points per 
event) for the "AND" events of React. [T] (139 events with TOF and <3>x cuts). The open 
histogram is the x p /xp distribution from inclusive single diffraction |2|, with a TOF veto 
only on the trigger side. The shaded and open histograms are normalized to the same 
area for the bin, 0.990 < x p < 0.995; (b) Same as (a), but the shaded histogram is x p or 
x p for the "OR" data (698 events with TOF and $x cuts). The vertical scale is linear. 
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Figure 15: Geometrical acceptance and track reconstruction efficiency vs. Mx- As de- 
scribed in the text, retention efficiencies of TOF veto cut and event selection are included; 
(a) "AND" triggered events withp andp both having 1.0 < — t < 2.0 GeV 2 . As discussed 
in the text, the open histogram assumes isotropic decay and the shaded histogram as- 
sumes longitudinal decay, for Mx > 4 GeV. The combination of solid and dashed curves 
is the acceptance function used in cross section calculations; (b) "OR" triggered events 
with 1.0 < — t < 2.0 GeV 2 for the observed final-state particle assuming isotropic decay. 
The solid line for Mx > 10 GeV is the assumed acceptance for longitudinal decay which, 
as for the "AND" data, is ~ 25% smaller than the acceptance for events with isotropic 
decay. 
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Figure 16: Monte-Carlo study of invariant mass calculation using the calorimeter. The 
observed (downward) shift in mass (true - observed) vs. the true mass. The fitted line 
corresponds to Eq. [8] in the text. 
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Figure 17: The mean "missing mass" calculated from the observed p and p, vs the 
corrected invariant mass calculated from the calorimeter information. See discussion in 
text. The vertical error bars on each point are the errors-in-the-mean for the missing 



mass calculation, while the horizontal bars show the event binning. As seen in Fig. |18 
there are zero events in the Mx range 14-18 GeV. 
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Figure 18: Final event sample; number of observed events vs. corrected calorimeter 
mass, M x , with 1.0 < -t < 2.0 GeV 2 ; (a) "AND" triggered data (85 events); (b) "OR" 
triggered data (586). 
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Figure 19: Mean number of struck calorimeter cells with at least 200 MeV energy (total 
energy in electromagnetic and hadronic sections) vs. corrected calorimeter invariant mass. 
Dashed line is the naive expectation, using < N >= 1.2M X , as discussed in the text. 
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Figure 20: Number of "AND" events with Mx < 6 GeV vs. ratio of electromagnetic 
energy to total energy detected in calorimeter. The peak corresponding to (e.m. energy 
= total energy) in both data and Monte Carlo is due to low-energy charged tracks losing 
all their energy in the e.m. sections of the UA2 calorimeter. 



43 




Figure 21: Central system decay distributions. dN/dcos(6) for all "struck" cells, averaged 
over the event sample: (a) for Mx < 5 GeV; (b) for Mj > 5 GeV. Histograms are the 
Monte-Carlo distributions described in the text. Vertical scale is arbitrary and linear. 
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Figure 22: daupE/dMx, the corrected differential cross section for React. [I] (proportional 
to the ratios of Figs. |TB| and Figs. |T3]), only for momentum transfer (s), t, of the observed 
trigger particle(s), p and/or p, in the range, 1.0 < —t < 2.0 GeV 2 ; (a) "AND" triggered 
data; (b) "OR" triggered data. The histograms are Monte-Carlo predictions assuming 
M^-independent, a^ v = 1 mb. As discussed in the text, the absolute values shown 
assume the (somewhat arbitrary) value, K = 0.74 GeV -2 . 
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Figure 23: Mass dependence of the Pomeron-Pomeron total cross section, derived 
from the "AND" and "OR" triggered data, respectively. The (arbitrary) cross section 
scale assumes K = 0.74 GeV -2 , as explained in the text. Dashed curve is the factorization 
prediction (which is independent of the assumed value of K) for the Pomeron-exchange 
component of Oy$p. The solid line is the fit to the "OR" points of a 7£eggeon-exchange 
term, (M^) -0 32 , added to this Pomeron-exchange term. 
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Figure 24: x p (xp) distribution for the "AND" data as in Fig. 0(a). The solid curve 
normalized to the data is the Monte-Carlo prediction assuming no explicit s'-dependence 
in c&L. 
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Figure 25: (a) Dashed curve is er-p*, the Pomeron-proton total cross section, from 
Ref. ||, assuming K = 0.74 GeV~ 2 (see text for explanation). The solid curve is only the 
Pomeron-exchange component; (b) same as (a), except for the proton-proton total cross 



section 34, 35 . 
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Figure 26: Predicted differential cross section daDPE/dM x (integrated over all t), assum- 
ing constant 0%%, = 1 mb for the Tevatron (y/s = 2 TeV) and the LHC (y/s = 14 TeV). 
In each case, the solid (dashed) curves are for assumed effective Pomeron intercepts, 
a(0) = 1.035 (1.00) respectively. 
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Figure 27: Predicted differential cross section d<jDPE/dM x (integrated over all t), for 
the HERA-B experiment with Pbeam = 920 GeV, assuming fixed-target pp interactions 
and constant a™ = 1 nib. 
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